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We study the formation of caustics in vortex-dominated flows. We find that only particles starting 
within a critical distance of a vortex which scales as the square roots of the particle inertia and the 
circulation can form sling caustics. We show that particles starting in an annular region around 
this critical radius contribute the densest clusters in the flow. The large density spikes occurring for 
such particles, even at small inertia, are indicative that these particles will experience large collision 
rates. 


I. INTRODUCTION 

When a turbulent flow contains small particles of a different density dispersed in it, these particles tend to cluster 
into small regions of the flow rather than being homogeneously distributed in space [lH5|. Several different but related 
ways of explaining this clustering have been posited Typically, heavy particles are centrifuged out of vortical 

regions and cluster into regions of high strain fTo|. This effect can observed for particles that are larger (typical of 
laboratory studies, see e.g. i) or smaller than the Kolmogorov scale. In particular, the behaviour of heavy particles 
much smaller than the Kolmogorov scale can be described using a simple equation (see eq. (1)) and is important in 
atmospheric flows. 

Particle inertia also leads to the formation of caustics. At a given time instant, the velocities of two neighboring 
particles can differ greatly from each other. Velocity contours of inertial particles may thus display discontinuities 
and multi-valuedness in some regions of the flow. Such regions are known as caustics, where particle velocity may 
not be described as a sing le-valued field. The effects of caustics have been extensively studied [llMl9l | . For example 
Mehlig and Wilkinson [ll| describe caustics as occuring due to folds in the velocity in phase space, leading to points 
in space where there are multiple values for the velocity at the same instant of time. Falkovich et al. [2Cj describe the 
same phenomenon as the formation of caustics due to the ‘sling effect’ or centrifugal action near the vortices. Caustics 
increase the collision rates for particles drastically, and have been suggested as a mechanism for the suddenness of 
the onset of rain in converting cumulus clouds [ll|, [lj, |20|• 

We ask how important the processes of sling caustics and preferential concentration are relative to each other. To 
this end we examine the archetypal flow for the sling effect: the flow due to a point vortex. In the process of being 
centrifuged out of vortical regions, particles which are closer to a vortex attain higher acceleration, and can eventually 
catch up with, or even overtake, particles which were initially further away from the vortex, thus forming caustics. 
We show that caustic-formation is only possible for particles initially located within a critical distance r cr ~ 1/2vTv 
of a point vortex, where T is the circulation of the vortex and r is the relaxation time of the particles, such that 
the particle Stokes number at the critical radius, St = Tr/r% r , is of order 1. Interestingly we find that the largest 
contributions to clustering come from the band of particles located in an annular region around this critical radius. 
This has consequences both for collision of particles and for condensation inhomogeneities. 

We note that the system of particles around a vortex has been studied previously by Raju & Meiburgj2l|. They 
provide a detailed discussion of the different particle-inertia regimes involved, and find similar spikes in density as we 
report in section Mil They also report an increased “accumulation rate” for intermediate values of the particle Stokes 
number. Their work was done before the idea of inertial particle caustics came about, and we find that caustics 
provide a way of understanding their results. 

In section an we discuss the general nature of the problem, including the different possible initial conditions that 
can be used. In section mu we look at point vortices in particular. We discuss when caustics can form in the flow 
around a vortex, and show some scaling laws that are obeyed. Since point vortices are idealisations and viscosity 
smears them out, we examine caustics formation due to Gaussian vortices in IIVI In section El we report simulations 
with many vortices in a doubly-periodic box. We conclude in section ED 


II. CAUSTICS AROUND A VORTEX 

We prescribe a two-dimensional steady axisymmetric flow field, with its centre at the origin of the coordinate 
system. The radius a of the particle is taken to be much smaller than the typical flow length scales. The motion of 
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such an inertial particle in a flow field u (x, t) is given by the Maxey-Riley equations [22|]. For particles much smaller 
than the Kolmogorov scale, and much denser than the fluid, the Maxey-Riley equations simplify to 


dr 



dv u — v 


dt 


(1) 


where r and v are the particle location and velocity respectively, and t is time. The effect of particle inertia is 
encapsulated in the particle time-lag r. 

The most general axisymmetric flow field is given by 
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where fl (r) is the angular velocity of a fluid particle at a radius r. For solid body rotation, p = 0, whereas for a flow 
with a single point vortex placed at the origin, p = 2. In cylindrical polar (r, 6) coordinates equation (jTJ) is written, 
with 


oj = dO/dt , 
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where the dot superscript stands for d(.)/dt. Particle d yna mics in the vicinity of one vortex has been studied using 
the above equations in some detail by Raju & Meiburg |2l| and by Shaw et al. [g|. 

From these solutions it is known that the radial distance of a particle from the center of the coordinate system is 
a monotonically increasing function of time. At long times the tangential velocity of a particle will approach that 
of the fluid, i.e., w ~ and so from equation ED, we have at long time or large radial distance, v r — rrfl 2 and the 
particle’s radial acceleration is negligible. 

We prescribe two rings of particles at initial radii roi and ro 2 with initial angular velocities LOi = coo(ri) , i = 1,2. 
Equations Q then give 


dAr 


d 

dt 


A 


dAv r 


dt 



dt 

Av r 

h- 

T 

(r 2 w) 


Av r 

A (ruj 2 ) 

A (r 2 fl) 
r 


( 5 ) 


where A(.) = (.)2 — (-)i represents a difference in a given quantity between the two rings of particles. In the following 
we define Aro = A r(t = 0) and use the shorthand notation ro for the initial location roi of the inner ring. 

A caustic can be said to have formed if, at a particular finite time t , particles on both rings arrive the same radius 
r but with different velocities, i.e. if Ar = 0, but Av r ^ 0. Note in particular that particles which asymptotically 
approach each other at long times do not give rise to caustics, although they do contribute to clustering. When a 
caustic occurs, trajectories in the phase portrait for the relative position and velocity look like the curve ‘b’ in figure 
|]]and not like curve ‘a’. 

We can solve the system ([5D for the general vortex S2 ~ r~ p and plot the separatrix in ro — Ar space between 
caustic-forming and no-caustic regimes. We only show results for p = 2, and comment on the implications of p ^ 2 
(see end of section ED. 

We use three kinds of initial conditions for this study. The first, one we shall call co-rotating (or ‘zero-inertia’), 
that the particles follow the flow exactly until t = 0, and particle inertia is switched on at t = 0. For particles of small 
inertia this is a reasonable condition |2 ll | . In the second initial condition, particles are initially motionless and the 
vortex is introduced at t = 0. We shall refer to this as the zero-velocity initial condition. A turbulent patch moving 
into a particle-rich region would correspond to this situation. Most studies of particles in turbulence employ one of 
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Figure 1: Schematic of possible trajectories for the system ©• At long times, the two rings of particles will come closer to 
each other and merge at infinite time, and this would, in the absence of caustics, give rise to a trajectory of type (a). In the 
caustics forming case, the ring closer to the origin at initial time would catch up with and overtake the ring further out, so the 
relative velocity would cross zero once in finite time, as in trajectory (b). 


these two conditions. Also consider a vortex tube being stretched by the turbulent flow around it. The strength of the 
vortex will grow with time until it reaches its maximum value. To mimic this situation, a third initial condition, which 
is just a generalisation of the second, is where both particles and vortex are initially motionless, and the strength of 
the vortex is gradually ramped up to its final value. The results from this third initial condition are not presented 
here because they display no unexpected physics. The introduction of another time scale of growing vorticity merely 
scales the answers appropriately. Most real-life situations will correspond with reasonable fidelity to one of our initial 
conditions or something in-between, so we expect our results to have a bearing on understanding clustering. 

The time and length scales inherent to the problem respectively are the inertial timescale r of the particle and 
L = vTV where T is the strength of the vortex. Note that T and r are the only parameters in the problem. For a 
given initial condition, from dimensional arguments we may argue that the caustic-formation time T caust i C will have 
to depend on nondimensional groups formed out of the scales of the problem, and must have a functional form 

F 'caustic rji ( r 0 Ar \ fn\ 

-^ = T \7rr’Wr)' (6) 

If there is indeed a regime where it is impossible for caustics to form, the caustic-formation time must diverge in this 
regime. We may then define a critical initial separation A r cr , as the maximum value of the initial separation for a 
given ro up to which caustics can occur. This quantity, suitably scaled, will be a function only of the initial distance 
from the centre of the vortex, i.e., 


A r cr _ ( r 0 \ 
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We study in detail the case of a single point vortex, and then examine the case of Gaussian vortices. Finally we 
will derive a necessary condition for caustics formation for a general axisymmetric flow. 


III. POINT VORTEX 


For a point vortex, we substitute p = 2 in equation @. The flow field is now given by 


u = 
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Figure 2: The separation between particle-pairs as a function of time, illustrating the difference between the types of initial 
conditions. The curves that do or do not cut the A r = 0 axis are respectively for particles that do or do not form caustics. 
Note that a maximum of one crossing of the Ar = 0 line can occur for a given curve. 


Using t and V^IV as scales, equations [4] may be written in non-dimensional form as 


x + x 
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where x = r/(rT) 1 / 2 and £ = r 2 w/F. We note that these equations are independent of r, i.e., the existence of caustics 
can be investigated with the particle size scaled out. From equation © it is clear that the azimuthal motion merely 
relaxes from any initial condition to (= 1 on a time-scale of r. At longer times, and when x 1, irrespective of the 
initial condition, equation ([5]) thus reduces to one where x is negligible and x ~ f{x). In such a regime, given that 
v r > 0 (because of the centrifugal nature of particle motion), it is obvious that two particles which are at different 
radii at a given time t in this regime will never produce a caustic. It follows that only particles which start out at 
small radial distance will accelerate enough to form caustics. 

A numerical solution of equation © enables us to divide the parameter space into regimes where caustics formation 
is possible and regimes where it is not. Figure [2] shows examples of both of these, and how the separation between 
the two rings of particles progresses with time. 


A. Limits for caustics formation 

The time at which each curve crosses the Ar = 0 line is T ca usUc- We expect the caustic-formation time to be small 
for small initial separations of the two rings of particles, and plot this quantity for A?’o = 10 -4 in figure [3] It is 
immediately clear that the caustic-formation time diverges even for this small separation above an Xq of about 0.5, 
for both initial conditions. We define the initial radius of the inner ring, above which caustics cannot form, as r cr . In 
the region where ro < r cr , we see that T caust i C scales as rg for the zero inertia initial condition, and as tq for the zero 
velocity initial condition. Note that, for the zero inertia initial condition, T caust i C is seemingly independent of r, since 
Tcaustic/ T scales as (ro) 2 /(Fr). After Raju & Meiburg|2l|. we can argue for this condition that the initial velocity 
term in equation [H] is negligible. This gives an equation for the evolution at early times of the separation between 
particles that goes as (assuming Ar <C ro): 




This suggests that T caU stics should scale as Tq for a given Ar. For the zero velocity initial condition, this simplification 
is not possible, and we can think of no simple argument at this point for the scaling as St -1 / 2 . 
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Figure 3: The caustic-formation time T caus tic as a function of the initial location of the inner ring of particles. A very small 
separation between the rings is used in this case, and the location of divergence is insensitive to further reduction in separation. 



Figure 4: Maximum initial separation for caustics formation for the zero-inertia and the zero-velocity initial conditions. Caustics 
can occur below of to the left of the lines. The curves for different r collapse onto each other when they are scaled according 
to equation [7) For the zero-inertia initial condition, ro and A r cr are inversely related. 


Figure |3] shows the maximum initial separation (Ar) cr for caustics to occur, as a function of ro- It is evident 
that beyond a certain ro = r cr , however small the initial separation, caustics cannot occur. For initial particle 
locations within a critical radius, one may have caustics, created with particles which are located within a certain 
critical separation distance (A r) cr . This quantity depends on the initial condition. For a range of ro, with initially 
co-rotating particles, the initial radius and the maximum separation for caustics obey an inverse relationship, i.e. 
ro(Ar) cr /(Fr) is constant. The zero-velocity initial condition, on the other hand, yields a much lower Ar cr which 
varies much more slowly with ro. 

If the inverse relationship of (A r) cr and ro is combined with the scaling argument for T caust i c , we can obtain a 
scaling for a "critical" velocity (A r) cr /T caust i C ~ TV/r^j. Comparing this with equation ([I]) suggests that the particles 
that bridge the largest gaps move outwards almost ballistically and the acceleration term is unimportant for such 
particles. 
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B. Particle density 

We know that inertial particles in a vortical flow cluster into regions of high density. We would like to understand 
the role of sling caustics in this process. We track the density of inertial particles using the Lagrangian particle 
tracking method of Osiptsov following Healy & Young j23j. This approach has been used in the past (see 0[25l). 
The equations are given in the appendix. By this approach, one may track the density in the vicinity of a given 
particle while moving with the particle. The method relies on keeping track of the net divergence of particle velocities 
in this moving vicinity. The Osiptsov method allows for obtaining local densities without actually tracking many rings 
of particles. Figure[5] (a) shows density profiles as functions of time in the vicinity of particles started at different ro. 
Densities in the vicinity of particles whose initial locations lie within r cr , i.e., those which can participate in caustics 
formation, are shown by dashed lines, whereas for particles starting out at ro > r cr , i.e., in regions which do not allow 
caustics, local densities are shown by solid lines. The difference between the two is stark. Outside the regime that 
allows caustics, all that happens is that particles cluster slowly with time, and the (Lagrangian) density monotonically 
increases. On the other hand, densities of particles starting from within the regime of caustics formation show a sharp 
spike around t = T caust i C for that ro, and subsequently the Lagrangian density decreases monotonically, since particles 
which participate in caustics formation later display trajectories that diverge in radius. The spikes happen by virtue 
of a patch of particles that start out within a small neighbourhood of each other clumping together. The most striking 
result in this figure is that the spike for ro ~ r cr is significantly taller than all the others, and the clustering remains 
higher for this case even at later times. We shall present further evidence below to show that the annular region close 
to the critical radius is the most important for clustering. 

In figure [5] (b), we show density profiles as a function of radius at a time of t = 15r and t = 20r. By this time, 
caustics have formed and the participating particles have since dispersed. In this and following figures, when we 
allude to "only caustics" or "excluding caustics", we refer to initial particle locations r 0 < 0.8r cr and r Q > 1.2r cr 
respectively. The values 0.8 and 1.2 have been chosen arbitrarily, and the answers do not change qualitatively when 
we choose other limits. We refer to initial locations within 0.8r cr < r 0 < 1.2 r cr as the annular region, which will be 
shown to contribute significantly to the densest regions. To identify where particles originated from, different symbols 
are given in this figure to the contributions of each of these. We see that there is a minimum radius below which there 
are no particles, i.e., the density is zero. This minimum radius increases slowly with time, at the rate of 1/a: 3 . We now 
discuss the large r portion of this figure. It is interesting to note that at a given radial location, two possible values 
of density are seen. This is because every radial location contains one set of particles which started at r* < r cr and 
another which started at r Q > r cr . Particles that start close to the vortex are thrown out so violently that they reach 
locations far away from the vortex. Their radial locations are now distributed over a much large radius. On the other 
hand, particles which started at r 0 3> r CT have moved much less and display a local density only slightly different 
from its initial value. Thus the neighbourhoods of caustics and non-caustics particles are of low and high density 
respectively. The picture changes at lower r. Here, the densest neighbourhoods are of the parcels which started out 
in the annular region near r cr . 


IV. GAUSSIAN VORTICES 

Viscosity acts on point vortices by smearing them out into Gaussian vortices. We therefore study how caustics form 
around Gaussian vortices. We do this by working with Gaussian vortices of fixed maximum vorticities and widths, 
such that their total strength is the same as for the point vortices. The simulations presented here are inviscid. 
Figures |B] plot the caustics boundary for the zero inertia and zero velocity initial conditions. We find that caustics are 
still created by particles initially within the region given by the same critical radius as with point vortices. However, 
particles which start out very close to the centre of the vortex do not participate in caustics formation in this case, so 
there is a minimum initial separation of vortex centre and particle for caustics formation. This radius is dependent 
on the Stokes number of the particles. This is because the effective value of the exponent p in eq. Ovaries with radius 
for a Gaussian vortex. Notice that the independence from particle inertia r of equations © and © is only valid for 
point vortices, i.e., when p = 2. For p values which are lower, we have made computations which show that the results 
are Stokes number dependent, and below a particular p for a given r, caustics cannot form. For any r, caustics form 
only over a range of ro between 1.2 r„ and r cr , where r v is the characteristic radius of a Gaussian vortex, defined by 

~ r / r v 


w(r) = oj v e 


( 10 ) 
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Figure 5: Lagrangian densities for initially co-rotating particles around a single point vortex at the origin. The densities 
were tracked using the method of Osiptsov (eq. 11411 . (a): Density as a function of time for different initial location tq. The 
caustics-producing initial conditions show sharp spikes, (b): Particle density as a function of radial distance from the vortex 
at t = 15r and t = 20r. There are two possible densities at each radial location, one in the vicinity of particles which started 
within the caustics-producing region and the other in the neighborhood of particles which started outside this region. 




Figure 6: Caustics with Gaussian vortices. These plots are analogous to the plots in figure [4] The maximum vorticity and 
characteristic radius are the same for all plotted curves. The dashed lines to the left of the curves show the characteristic 
radius of the Gaussian vortex normalised by the largest inertia. Caustics exist only for r cr > ro > 1.2 r v , and the inner limit 
is dependent on St (see text for a definition of St). Left: zero-inertia initial condition, right: zero-velocity. The values of St 
shown are 0.04, 0.08, 0.12, 0.16, 0.2 and 0.4. 


uj v being the peak vorticity. The angular velocity of a particle around the origin is therefore given by 

fi(r) = J 2nrudr = (l - e~ r2/r = | (l - e~ r2/r . (11) 

The T (= 7roj v Ty) in equation fill is strength of the Gaussian vortex. Figures [6] show that the vortex strength plays 
much the same part for Gaussian vortices as it does for point vortices. The Stokes number of a particle may then be 
















































defined using the characteristic radius of the vortex: 


St = - (12) 

1 V 

Figures [G] plot the analogues of figure 01 The actual value of r v does not make a difference to critical radius r cr . 
However, the nondimensional parameter lo v t (which is a multiple of the particle Stokes number and hence depends 
on r v ) limits the occurrence of caustics. No caustics occur for lo v t < 4. Therefore, as the width r v of a Gaussian 
vortex increases in time due to viscosity, only particles of larger inertia can undergo caustics at later times. 


V. SEVERAL VORTICES 


We now apply the same density tracking procedure as in sections 1III1 and IIVI to a system of particles in the flow of 
several point vortices in a periodic box. The system of point vortices is for us a proxy for two-dimensional turbulence. 
(While Gaussian vortices would have been more realistic, no expression exists of the form [T3] for Gaussian vortices in 
a periodic box.) 

We perform simulations with N = 10 point vortices of the same circulation but random signs (zero net circulation) 
in a doubly periodic box. The motion of the vortices is computed using the Hamiltonian for this problem [26| : 



^(Im {Zj k )f , 


(13) 


where Zj is the position of the jth vortex in the complex plane: Zj = Xj + iyj ; L x is the length of the box along the 
x-direction and A = L x L y is the area of the box. 

We study statistics of the particle density as a way of quantifying the effects of caustics in such systems. We 
perform simulations of three types: (i) by distributing particles at random but removing all particles inside the caustic- 
producing region of each vortex. In simulations of type (ii) we do the reverse, i.e., only keep particles within the critical 
r cr for caustics formation. In simulations of type (iii), we start particles in the annular region 0.8r cr < r < 1.2r cr . 
Simulations excluding caustics are done with 4 x 10 4 particles; and simulations with particles only in the caustic 
regions, and particles in an annulus are done with 10 4 particles each. We hold r constant during each simulation. 
Figure [3 shows particle density in the domain at different times, for initial particle locations exclusively within, and 
strictly outside, the caustics-forming region, i.e., for particle distributions corresponding to (i) and (ii) above. We 
notice that the highest densities occur in both cases in the same regions. These appear to correspond to initial 
conditions at the edge of the caustics-forming region, i.e., around the critical radius. To confirm this, we perform 
simulations of type (iii). Figures [8] show how significant the region around the critical radius is for producing high 
density particles, in comparison to the regions well inside and well outside. We also note that the effects of the critical 
annulus region are more pronounced for the smaller of the two inertia values. 

Another point of interest is that there appear to be, in figure [3 "voids" in the distribution of particles in space, 
with the radius of the voids given by the critical radius. It is known from simulations of two- and three-dimensional 
turbulence ana that such voids exist. It is also known that, for small inertia, the voids increase in size with 
increasing inertia. It is subject to future confirmation whether our idea of a critical radius r cr for inertial particle 
caustics is directly connected to these voids (and the dependence of their size on particles inertia). The fact that 
caustics formation is much faster than flow time-scales, and that a given particle inertia leads to a characteristic void 
size lends credence to this. 

The number of particles with greater than a threshold density, as figure [3 shows (we have checked that the exact 
value of the threshold does not matter; for this figure a value of 1000 was used) is also much higher for particles 
starting in the annulus. In fact for the smaller of the two inertias, almost all ‘caustic events’ happen for particles 
starting in the annulus. 

We plot the distribution of Lagrangian densities for the many-vortex simulations. The plots for particles starting 
in the ’only caustics’ and ’excluding caustics’ corroborate the results from section EU figure [3 particles starting 
outside the caustics region show a monotonically increasing density in their vicinity, whereas particles starting close 
to vortices get thrown out violently and so a peak in the distribution is seen at very small density. In figure 1101 we 
plot the density distributions for particles starting in the annulus regions for two different values of particle inertia. 
The curves are plotted at the same advection time (for identical fluid flows, since the particles don’t affect the flow), 
which is a different multiple of each different inertia. The curves show the smaller of the two inertias reaching a higher 
average density at the same advection time, which is expected, since the heavier particles have had less time in units 
of t to cluster. 
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Figure 7: Snapshots at time t = 200r of particle positions with 10 vortices in a doubly periodic 10 x 10 box. The colour 
indicates local particle density. The particles have t = 0.001. Left: Caustics excluded; right: Only caustics. 
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Figure 8: Average Lagrangian densities. The three curves in each plot correspond to particles starting well within, well outside, 
and in the annulus around the critical radius r CT . Particles starting close to r cr are the primary contributors of high density 
regions. Left: r = 0.001, Right: r = 0.01 


VI. CONCLUSION 

To summarise, we have studied the formation of caustics in model two-dimensional vortical flows. For a single 
vortex, we show that only particles which are initially within a critical radius r cr ~ 0.5( Ft ) 1 / 2 from its centre can 
cause caustics to form. For a point vortex, particles placed inside this critical radius can only overtake rings of 
neighbouring particle rings which lie within some initial separation (A r) cr . The time for such overtaking is the 
caustics formation time T caust i CS - This quantity diverges for ro —> r cr . The initial separation within which caustics 
form depend on initial conditions and so does the formation time. 

The formation of caustics is intimately related to the creation of regions of high particle density. Our adoption of 
Osiptsov’s method enables us to track densities in the neighbourhood of Lagrangian particles . Interestingly, particles 
which lie close to the critical radius, i.e. on the edge of the caustics-producing region, give rise to the largest spikes 
in density, suggesting that this region should contribute to the effects that particle caustics have on the flow. Indeed, 
we find that for small inertia, almost all the instances of very dense clusters come from the annular region. In these 
regions particle collision frequency is likely to be the highest . Further studies are needed to confirm this prediction. 
Outside the critical radius, particle density displays a qualitatively different behaviour. Density increase is gentle and 
monotonic in this case for all time. 

Our results above show that the centrifugal force can only create caustics in a small area around a vortex. We 
also show that the biggest density spikes happen not for particles that start well within the caustic regions, but rather 
for particles that start on the edge of the caustic regions (i.e. for particles which may be considered to have Stokes 
numbers of order 1. We find also that the effects of this annular region around the critical radius are more prominent 
for the smaller of the two inertias. 

Our results predict that particle collisions in flows with strong vortices are more probable at a critical radius from 
the vortex. Our future work will test this prediction with more realistic flows than a system point vortices. Whether 
the existence of this annular band depends on the confinement offered by the two-dimensional nature of our setup 
is also worth checking. We hope that our work will motivate viscous three-dimensional simulations to check our 
predictions on the importance of the critical annular region about each vortex in the creation of the densest clusters. 
In particular it would be useful to ask whether the extent of clustering can be directly correlated to suitably weighted 
summations of square roots of instantaneous circulations in the flow. 
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Figure 9: Fraction of particles in high-density regions. We plot the fraction of particles whose neighbourhoods have density 
greater than a certain threshold, here prescribed to be p cu t 0 ff = lOOOpo, where po is the uniform initial density. Left: r = 0.001, 
Right: t = 0.01. 



Figure 10: The density distribution for particles starting in the annulus regions for two different values for inertia. The curves 
for the larger inertia seem to shift towards lower density values. 
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Appendix: Lagrangian Density tracking 
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The density tracking equations of Healy & Young 
follow the same notation as them. They are 


are modified for a cylindrical polar coordinate system. We 
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The equations |T7] can be integrated in time to find the density given by 
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(15) 


The first of equations [15] includes the factor r because of the Jacobian involved in transforming from Cartesian to 
polar coordinates. 

A comment on the validity of Lagrangian density tracking Osiptsov’s method as described above calculates 
the density around a particle as it moves along on its trajectory, and does this without reference to how many other 
particles are around the particle being tracked. We have found that Lagrangian density tracking is better than the 
usual way of calculating density (which is to count the number of particles per unit area/volume). However, the 
density predicted by the method will only be realised in a real flow if there is in fact a continuum of particles that 
converge onto (or diverge away from) the particle being tracked. Since this is hardly observed in real flows (in clouds, 
for example, there is on average one particle per every few Kolmogorov boxes), the answers from Osiptsov’s method 
must be taken with a pinch of salt. 
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